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EVAPORATION,  HEAT  TRANSFER,  AND  VELOCITY  DISTRIBUTION 
IN  TWO-DIMENSIONAL  AND  ROTATIONALLY  SYMMETRICAL 
LAMINAR  BOUNDARY-LAYER  FLOW* 

By  Nils  Frossling 


INTRODUCTION 


Aside  from  the  simple  case  of  the  plane,  no  quantitative  calcula¬ 
tions  of  the  evaporation  of  a  body  in  a  moving  medium  exist  so  far.  The 
heat  transfer,  which  under  certain  circumstances  (see  p.  4)  follows  the 
same  laws,  has  been  treated  theoretically  for  the  cylinder  by  Kroujiline 
(ref.  1)  and  Squire  (ref.  2).  For  boundary-layer  flow,  Kroujiline  used 
for  the  temperature  field  a  power-series  method  of  the  same  type  which 
has  been  introduced  for  velocities  by  Pohlhausen  (ref.  3).  Because  of 
this  stipulation  of  the  profile  form,  the  result  must  be  approximate, 
and  the  eventual  agreement  with  the  correct  value  is  rather  accidental. 
Squire  gave  an  exact  treatment  of  the  transfer  in  the  immediate  proximity 
of  the  stagnation  point.  However,  it  is  of  great  interest  to  have  a 
calculation  method  for  heat  and  mass  transfer  in  the  entire  boundary 
layer,  the  error  of  which  depends  only  on  the  work  expenditure  of  the 
numerical  calculation  and,  therefore,  not  on  possible  approximative 
formulations.  Even  though  the  calculation  is  time  consuming,  one  has 
the  advantage  of  being  able  to  check  approximate  and  more  rapid  methods 
with  respect  to  this  solution.  Under  the  supposition  that  the  constants 
of  the  problem  (shape  of  body,  pressure  distribution,  etc.)  may  possibly 
be  eliminated  from  the  equations  to  be  solved,  it  is  also  possible  in 
several  cases  to  vise  the  complete  exact  solution  directly.  The  author 
of  this  report  perfected,  for  this  reason,  the  exact  solutions  for  the 
temperature  and  concentration  fields.  Two  dimensional  and  rotationally 
symmetrical  steady  boundary -layer  flows  were  treated.  The  latter  case 
is  the  more  complicated  one  because  of  the  form  of  the  continuity 
equation. 


*"Verdunstung,  Warmeubergang  und  Geschwindigkeitsverteilung  bei 
zweidimensionaler  und  rotationssymmetrischer  laminarer  Grenzschicht- 
strdmung."  Lunds  Universitets  Arsskrift,  N.F.  Avd.  2,  Bd.  36,  Nr.  4 
Kungl.  Fysiografiska  Sallskapets  Handlingar,  N.F.  Bd.  51,  Nr.  4,  1940. 
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For  the  calculation  of  the  transfers,  however,  the  velocity  fields 
must  be  known.  For  the  two-dimensional  case  one  has,  aside  from  approx¬ 
imate  methods  of  solution  by  Pohlhausen  (ref.  3),  K&m&n  and  Millikan 
(ref.  4),  and  others,  the  method  of  Blasius  (ref.  5)  and  Hiemenz  (ref.  6) 
which  was  improved  by  Howarth  (ref.  7)*  By  means  of  this  method  one  may 
solve  the  equations,  without  any  arbitrary  assumptions  regarding  velocity 
profile  and  the  like,  by  power-series  development  from  the  stagnation 
point  up  to  an  arbitrary  point  on  the  meridian  curve.  The  work  expend¬ 
iture  depends  on  the  required  accuracy  and  on  :,ne  position  of  that  point. 
Since  the  development  becomes  very  rapidly  more  cumbersome  with  the 
distance  from  the  stagnation  point,  it  is  appropriate  to  use,  from  a 
certain  point  onward,  a  continuation  method,  for  instance,  according  to 
Prandtl  (ref.  8)  and  Gortler  (ref.  9).  Howarth  (ref.  J)  transformed 
the  functions  of  Blasius  and  Hiemenz  into  functions  of  such  a  type  that 
the  constants  disappear  from  the  equations  so  that  the  numerical  solu¬ 
tions  for  them  may  be  applied  to  any  two-dimensional  flow.  He  treated 
the  symmetrical  as  well  as  the  unsymmetrical  case,  and  indicates  the 
solution  of  one  of  these  functions.  For  the  present  investigation, 
Howarth' s  functions  are  used  in  the  two-dimensional  case.  Since  the 
accuracy  of  Howarth' s  numerical  tables  is  not  sufficient  for  the  cal¬ 
culation  of  the  transfer,  a  new  numerical  calculation  was  made  of  certain 
functions.  For  the  three-dimensional  rotationally  symmetrical  case,  in 
contrast,  there  exists,  so  far,  no  calculation  of  the  functions  of  the 
series  development.  Because  of  the  modified  continuity  equation,  other 
systems  of  equations  must  be  used,  and  these  systems  are  established 
here.  Hie  necessary  functions  are  numerically  calculated.  Although  the 
form  of  the  meridian  curve  takes  effect  through  the  continuity  equation, 
one  can  proceed  in  the  distribution  of  the  functions  in  such  a  manner 
that  the  constants  of  the  meridian  equation  disappear,  and  the  solutions 
therefore  are  valid  not  only  for  arbitrary  pressure  distribution  but  also 
for  arbitrary  shape  of  the  body  of  revolution.  The  continuation  method, 
beginning  with  the  limit  of  validity  of  the  broken-off  power  series,  has 
been  perfected  also  for  this  case. 


An  investigation  is  carried  out  regarding  the  validity  of  the  law, 
stated,  for  instance,  by  Ulsamer  (ref.  10)  that  the  Nusselt  number  is 
proportional  to  the  cube  root  of  the  Prandtl  number. 

A  few  approximation  methods  for  the  calculation  of  the  transfer 
layer  are  discussed. 

Only  a  brief  survey  is  presented  here  since  a  more  detailed  report 
is  to  be  given  in  a  later  paper.  - 
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THE  FUNDAMENTAL  EQUATIONS 


With  dimensioned,  quantities,  the  boundary -layer  equations  for  flow, 
concentration,  and  temperature  read 


du  ,  du 
dx  dy 


UU'  +  v 


d2u 

dy2 


du  |  dv 
bx  by 


0  two  dimensional 


I  d(ur)  d(vr) 

<  - +  -  =  0  rotationally  symmetrical 

I  bx  by 


u 


be 

bx  V 


U  +  V 
dx 


dc 

dy 

dy 


A 


ay2 


a 


£i. 

9y2 


Concerning  the  derivation  of  the  two  last  equations  for  the 
rotationally  symmetrical  case  see  p.  15-16.  The  boundary  conditions  are 


y  =  0;  u  =  v  =  Oj  t  =  t0i  c  »  cmj 
y  =  oo;  u  =  U;  t  =  Oj  c  -  0. 


In  these  equations  the  customary  designations  of  the  various 
quantities  are  used,  x  =  distance  along  the  body  surface  from  the 
stagnation  point  to  the  base  point  of  the  normal  to  the  body  surface, 
y  =  length  of  that  normal,  r  =  distance  from  base  point  to  axis  of 
rotation,  u,  v  =  velocity  components  in  the  direction  of  x  and  y, 
respectively.  U  =  the  velocity  component  parallel  to  the  body  surface 
immediately  outside  of  the  boundary  layer  (calculated  from  the  experi¬ 
mentally  determined  pressure  distribution),  t  =  excess  of  the  temper¬ 
ature  of  the  surface  over  the  temperature  of  the  undisturbed  fluid, 
c  *  corresponding  concentration  quantity,  v  =  kinematic  viscosity, 

a  =  temperature  diffusivity 


A  =  diffusion  coefficient.  If 
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Uq  is  the  undisturbed  velocity  and  D  a  characteristic  length  of  the 

body  (for  instance,  its  diameter),  one  can  transform  the  equation  into 
dimensionless  form  by  dividing  the  velocities,  lengths,  temperatures, 
and  concentrations  by  the  quantities  Uo,  D,  to,  and  cm.  These  equa¬ 
tions  contain  as  constants,  among  others,  the  Reynolds  number  Re  =  — . 

v 

This  number  one  may  eliminate  by  modifying  the  scale  of  the  boundary 
layer  in  transverse  direction,  by  multiplying  the  values  of  y  and  v 

by  \JlRe.  The  equations  used  below  with  the  designation  "dimensionless 
equations  without  Reynolds  numbers"  are  changed  In  their  appearance, 
compared  to  those  mentioned  above,  only  in  that  v  disappears,  and  A 

and  a  are  replaced  by  &  and  In  the  boundary  conditions  to  and 

/\  £ 

cm  are  replaced  by  1.  The  two  quantities  —  and  —  which  are  often 

independent  of  pressure  and  temperature,  as  in  the  case  of  ideal  gases, 
are  dependent  on  the  media  used.  These  quantities  are  called  Stanton's 
numbers.  Frequently  their  inverse  values  are  used,  designated  as  Prandtl 
numbers.  In  an  earlier  report  of  the  author  on  the  evaporation  of  drops 
(ref.  11)  the  designation  a  is  used  for  the  Stanton  number.  Since  at 
present  this  letter  is  used  mostly  for  the  Prandtl  numbers,  this  defini¬ 
tion  is  employed  in  the  present  report  to  prevent  misunderstandings. 

Thus  a  here  signifies:  o  =  ^  or,  respectively,  =  ~  =  V<^>P» 

The  equations  for  temperature  and  concentration  are  therefore 
identical  when  t  and  c  are  interchanged.  For  the  temperature¬ 
boundary  layer  it  is  assumed,  however,  that  the  dissipation  and  the 
heat  generated  by  change  in  pressure  may  be  neglected.  This  assumption 
is  satisfied  for  not-too-large  velocities  (ref.  12).  The  equations  also 
presuppose  that  the  velocities  be  small  compared  to  sonic  velocity  in 
order  to  make  the  compressibility  negligible.  A  further  limitation  of 
the  equations  is  given  by  the  fact  that  the  differences  in  concentration 
and  temperature  must  not  be  so  large  that  the  constant  characteristics 
of  the  media  vary  from  point  to  point.  Because  of  the  identical  form 
of  the  two  equations  for  temperature  and  concentration,  which  is  thus 
satisfied  under  these  presuppositions,  both  may  be  treated  simultaneously. 
In  the  following  equations  one  may,  therefore,  immediately  interchange 
the  quantities  c  and  t. 

In  the  search  for  a  solution  which  satisfies  the  accuracy  require¬ 
ments  discussed  in  the  introduction,  the  method  of  power-series  develop¬ 
ment  in  x  was  used.  Breaking  off  the  power  series  after  a  certain 
number  of  terms  one  was  able  to  use  this  solution  from  the  stagnation 
point  up  to  a  point  the  position  of  which  was  dependent  on  the  accuracy 
requirements.  Starting  from  this  point  one  could  then  use  for  the  layers 
of  different  types  step-by-step  continuation  methods. 
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For  the  sake  of  brevity,  we  shall  use  below  the  common  name 
"transfer"  boundary  layer  for  the  temperature  and  for  the  concentration¬ 
boundary  layer. 


POWER-SERIES  DEVELOPMENTS  IN  x 


A.  TWO  DIMENSIONAL  CASE 
a.  Flow  Boundary  Layer 
1.  Symmetrical  case 


As  was  mentioned  in  the  introduction,  this  case  has  been  treated 
by  Blasius,  Hiemenz,  and  Howarth.  This  report  uses  for  the  most  part 
the  same  designations  as  Howarth.  The  only  difference  is  that  Howarth' s 
quantities  Fv  are  here  replaced  by  the  quantities  \|ry  because  the 

capital  letters  are  more  suitable  for  the  functions  of  the  transfer 
boundary  layer. 

In  order  to  replace  the  two  unknown  quantities  u  and  v  by  a 
single  one  (ijr),  the  following  conditions  satisfying  the  continuity 
equation  are  set  up  as  usual: 


u  = 


d\|r 

$y 


c)\|r 

dx 


The  first  flow  equation  then  becomes  (with  dimensionless  quantities 
without  Reynolds  numbers) 


&  _  b^  _  yy,  + 

by  bxby  bx  $y2  $y3 

Blasius  and  Hiemenz  solved  this  equation  by  means  of  the  formula 


for  the  symmetrical  case  where  the  velocity  distribution  outside  of  the 
boundary  layer  follows  the  formula  U  =  u^x  +  u^x^  +  u^x^  +  .  .  .  .  Vv 
are  functions  only  of  y. 
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Ejy  comparison  of  the  various  powers  of  x  the  equations  for  \|fv 
were  obtained.  These  equations  were  freed  of  the  constants  Uy  by 
introduction  of  the  functions  fv,  8v>  etc.,  by  means  of  the  following 
statements: 


2.  Unsymmetrical  case 

For  an  unsymmetrical  two-dimensional  body  for  which  the  velocity 
distribution  follows  the  formula  U  =  u^x  +  Ugx^  +  u^x^  +  .  .  .  ,  the 
formulation 


♦ 


♦ix  +  ^2x2 


+  •  •  • 


was  used.  Here  also  the  i/v  were  freed  of  the  constants  Uy,  this  time 
by  the  expressions 


n 


etc. 


b.  Transfer  Boundary  Layer 
1.  Symmetrical  case 

Uhe  author  of  this  report  attempted  in  the  dimensionless  equations 
without  Reynolds  numbers,  aside  from  the  formulations  for  \|r  mentioned 
above,  a  development  for  c  in  the  following  manner  (cv  are  functions 

of  y  only): 


£ 

v=0 


CyXV  = 


'0 


+  CiX  + 


C2X£ 


+  QxX' 


c  * 


+ 


7 


Boundary  conditions: 


,v. 

•  •• 

•• 

: 

y 

=  0 

c0  =  1 

o 

ii 

• 

• 

• 

it 

CM 

O 

II 

r H 
O 

:•••%  < 
•  • 

y 

=  00 

c0  =  C1 

=  =  •  •  •  —  0 

Bty  substitution  one  obtains  for  the  and  c2n+i 


le¬ 


an 


=  5Z2lD|f'2n+l-2kc2k  -  ^  (2n  +  1  -  2kH2n+l-2kc  '2k 
k=0  0 


1  c" 
a  c  2n+l 


n 


n 


-  £I(2k  +  W2kc2k+l“SI(2n  +  1  -  2k)* 


k=0 


0 


2n+l-2kC ' 2k+l 


One  can  easily  show  that  the  equations  for  are  such  that 

they  become  identically  zero.  In  the  groups  of  equations  mentioned 
above  which  constitute  the  recursion  equations,  there  occur  exclusively 
functions  with  even  or  odd  subscripts.  Prom  this  one  can  see  that  c 
Is  an  even  function  of  x  which  follows,  besides,  from  the  nature  of 
the  problem.  In  order  to  be  free  of  the  constants  Uy,  new  functions 

are  introduced  by  the  following  statements: 


_  ^uj5F2 

n  =  y  }  *v  as  above>  co  =  Fo»  c2  -  > 


\  +  N 


6  u-l  h  UluT  H6  Ui2^  61 


Uj 


10u9  L  .  u3“7  „  .  V  „  +  VS  J  + 

hQ  +  ~  Jg  +  ~ 


c8 3  ~^rv& 


Gq  +  — —  —  Hq  +  ■  ■  —  IYq  t  — 

1  \  8  uiu9  8  ulU9  8  ux2u9 


uiu9 


% 


Boundary  conditions: 

q  =  0;  Fq  =  1;  the  remaining  functions  =  0 


t|  =  “>; 


all  functions  =  0 
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The  remaining  equations  do  not  have  explicit  solutions  and  must  therefore 
be  integrated  by  other  methods,  for  instance,  according  to  Runge  and 
Kutta.  (See  p.  19.) 


2.  Unsymmetrical  case 

With  c  =  cQ  +  c^x  +  c^2  +  c^x3  +  .  .  .  and  the  boundary  conditions 


y  “  o;  co  *  1?  ci  -  c2  =  .  .  .  =  0 


y  =  ooj  c0 


—  C i  —  Co  —  •  •  •  0 


one  obtains  for  c. 


±cv"  = 

0  v-n 


c  - 


fco  v+1“k  k 


f>  +  1  '  l)Wkck' 


Here  none  of  the  functions  cy  disappears,  and  c  is  therefore, 

as  had  been  expected,  an  even  function.  The  equations  here  are  not 
divided  into  two  independent  groups  but  the  functions  follow  successively 
one  from  the  other. 


Distribution  of  the  cv: 


3up 

n  =  *1 =  fi\/^  etc,?  c°  =  F°;  ci  =  lJ 


4us 


C2  =  ^  c3  ■  Sf 


\  e,  .  ™tL  ,  1*2*  + 


ul% 


^2 

+  -== — K3 
u12u4 


6uc  /  UoUh  Ui^  up2u,  u2*t  \ 

— ■ 11  (hi.  +  -Z—Bl  +  — — +  -4—^4  » 

V5  ulu5  U1  u5  u15u5  / 


C),  = 


u 
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Boundary  conditions: 


t)  =  0;  Fq  =  1}  the  remaining  functions  =  Oj 
t}  =  °o;  all  functions  =  0; 


w 


If  " 

a  x 

IGo" 


1h2" 


X" 

W 


lv-  " 

X” 

X" 


ij4" 
a  H 


X” 


-W  +  fl’Fl  -  2*2*0 ' 

-f^’  +  2f1’G2  -  3g3F0’ 

► 

-flV  +  2VH2  -  »}V  +  f (f 2 ' Fl  - 

-fl°3'  +  3fi’<b  -  HFo' 

-flH3’  +  Siyilj  -  4h4F0 '  +  ^(2f2'G2  +  gj’F!  -  2f2G2'  -  3gjF1  'j  . 
-fiV  +  Jfi’Kj  -  AV  +  fl21^2  +  hj-Fj.  -  SfgV  -  31^ 

-fiV  +  4fiA  -  AFo' 

_flV  +  4flA  -  Vo'  +  !(5f2 'G3  +  <VF1  -  2f 2°3 '  -  Vi) 

-fA'  +  4V\  -  5Vo'  +  l(283,02  -  3g3a2') 

-f  x  ’  +  ,4f1'J^  -  5J5F0'  +  2^f2E3  +  VF1  '  2f2H3'  "  ^h4Fi  ^ 

y^g5’H2  +  h3'G2j  -  8^g3H2-  +  h3G2') 

-f^V  +  \  -  5q^FQ'  +  §(3f2'K3  +  k)+'F1  -  2f2K3'  -  kls^F^  + 

|(a3'H2  -  3h5H2i) 


The  first  equation  is  identical  with  the  first  one  of  the  symmetrical 


case. 
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B.  ROTATIONALLY  SYMMETRIC  CASE 
a.  Flow  Boundary  Layer 


For  flows  about  a  blunt  body  of  revolution  whose  axis  lies  in  the 
direction  of  the  flow,  the  flow  equations  are,  according  to  Boltze 
(ref.  13), 


du 
u — 

dx 


+ 


+ 


d(ur)  +  d(vr)  =  Q 
bx  by 


In  the  transformation  to  the  dimensionless  form  without  Reynolds 
number  v  disappears.  The  quantity  r  then  must  have,  for  the  bodies 
of  revolution,  the  meaning,  distance  of  axis  of  rotation  up  to  the  base 
point  of  the  normal  instead  of  up  to  the  point  (x,  y) .  A  function  for 
identical  solution  of  the  continuity  equation  is  desired.  Boltze 
(ref.  13)  suggested  a  function  \|f  which  is  defined  as  follows : 

1  dty  _  1  djr 

r  dy  r  dx 

In  the  present  report  another  solution  \(r  also  has  been  examined. 
Definition 

cty  1  b(^r)  cty  i|f  dr  1  d(tr) 

by  r  by  bx  r  dx  r  bx 


The  function  has  the  advantage  that  the  equations  become  some¬ 
what  simpler  and  that  the  velocity  profile  is  obtained  directly.  The 
functions  are  derived  in  both  cases .  and  i|rr  are  flow  functions . 


1.  Use  of  the  function  \|r 

After  substitution  of  the  expressions  for  u  and  v  into  the 
first  boundary-layer  equation  one  obtains 

1  +  r  §£  _  r  =  Puu*  +  r2 

\pyj  dx  by  bxby  bx  'dy3 
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The  power-series  developments  used  function  of  y;  rv  and 

uv  constants ) 

\Jf  =  +  tygx^  +  .  .  • 

3  5 

r  =  r-jX  +  +  r^x-'  +  .  .  . 

U  =  U-jX  +  u^x^  +  u^x^  +  .  .  . 


The  functions  4ry  have  the  following  boundary  conditions  . 


y  =  o;  -  V  =  °; 

y  =  ^2  =  rlul>  ty'  =  riu3  +  r3ui»  *6'  =  riu5  +  r3u3  +  r5uli  • 


After  substitution  of  the  power  expressions  into  the  equation  for 
\|r  one  obtains  equations  for  \(rv  by  comparison  of  the  different  coef¬ 
ficients.  These  equations  may  be  freed  of  the  letters  rv  and  Uy 
by  the  following  formulations: 


The  new  functions  have  the  boundary  conditions 

q  =  0;  all  functions  and  their  first  derivative  =  0 

T]  =  °°J  ?2  =  x5  64 '  =  h4'  =  66'  =  h6'  =  J6'  =  y  V  =  %'  =  °J 


0 


One  obtains  the  following  equations: 


«•  •• 

•  •  • 

•  •• 

•• 

•  • 

•  ••• 

•t  •• 

•  •  • 

•  • 


fa1"  -  -f2f2"  +  |(f2'2  -  1) 

64’"  =  ^264”  +  2VS4'  -  2f2,'B4  -  1 

V"  =  -f2h4"  +  2t2\'  -  2f2'\  -  4(5f2'2  -  2f2f2”  +  *) 

■H"  -  -f286"  +  3f2'86'  -  3f2"g6  -  1 

V"  -  -W  +  5f2'h6'  -  5f2"l>6  -  | (5f a ■ 2  -  Sfgfg"  +  1) 

kg'"  .  -f2k6"  +  3f2'k6'  -  3f2"k6  +  f(384 '2  -  H«U")  '  2 

V"  -  -faJ6"  +  3W  -  5f2"J6  +  H  V  -  f  (g4h4”  4  W') 

3  2g4  “  ‘‘W  +  2f2  g4  ■  ■*•) 

V"  =  -f2qg"  +  3f2'«Jg'  -  3f2"qg  +  -  to^")  +  §(fgh4 

^2%  +  ^(5f2'2  *  2f2f2") 


2.  Use  of  the  function  i|r 


One  obtains 


d\jf 

by  bxby 


d\|r  d^\|r 

bx  ^y2 


1  dr  —  ^2^ 

- \|r  — 5- 

r  dx 


=  UU* 


+ 


by3 


Power  series  developments 

i  =  tjX  +  +  \jf^X5  +  .  .  . 

3  5 

r  =  r-jx  +  r5x^  +  r^xy  +  .  .  . 
U  =  u-jX  +  u^x^  +  u^x5  +  .  .  . 
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In  the  first  method,  the  functions  have  even  subscripts,  in  the 
second  odd  ones.  A  simple  relation  exists  between  the  two  groups  of 
functions  which  one  may  easily  obtain  by  equating  the  two  expressions 
defining  u_  and,  respectively,  the  two  expressions  defining  v  (which 
gives  if  =  yr). 


f2  =  fi; 

fl 

g4  _  &y  -  ^5  +  g  * 

f 

g6  =  g^j  h6  =  h5  +  kg  =  k5;  J6  =  +  |  g3J  qg  =  +  |  hj 


b.  Transfer  Boundary  Layer 

Hie  general  equation  of  the  temperature  and  concentration  fields 
for  rotationally  symmetrical  flow  has  not  been  set  up  before.  For  the 
special  case  where  the  body  is  a  sphere,  the  author  (ref.  11)  has  shown 
that  for  boundary -layer  flow  the  equation  is  identical  with  the  one  for 
the  two-dimensional  case,  at  least  for  points  which  do  not  lie  directly 
at  the  stagnation  point.  In  the  present  report,  it  is  shown  that  the 
same  boundary -layer  equation  is  valid  also  for  arbitrary  blunt  bodies 
of  revolution,  and  that  this  applies  to  points  directly  at  the  stagnation 
point  as  well.  The  introduction  of  mass  or  heat  into  a  volume  element 
by  diffusion  and  convection  is  expressed  by  the  following  equation  (which 
is  valid  for  rotationally  symmetrical  flow  without  neglect  of  the  boundary 
layer  when  x  and  r  are  counted  up  to  the  element  instead  of  to  the 
base  point) : 


d(cur)  d(cvr)  fb  x,  \  fx>  \ 

-5T+-TTm*?:S7 


The  derivation  becomes  the  simplest  if  one  chooses  as  the  volume  element 
an  element  bounded  by  two  meridian  planes,  two  surfaces  x  =  constant, 
and  two  surfaces  y  =  constant.  In  order  to  arrive  at  the  boundary -layer 
equation,  one  groups  the  derivatives 

+  dc  1  dr  +  d2c  +  1  dr  dc 
dy^  dy  r  dy  dx2  r  dx  dx 


S.  d(ur),  +  c  d(vr )  +  udc  +  ydc=A 
r  dx  r  dy  dx  dy 


1 6 


Supposing  that  a  thin  boundary  layer  exists,  the  terms  2  and  3  in 
the  parenthesis  disappear.  The  first  two  terms  of  the  left  side  disappear 
because  of  the  appearance  of  the  continuity  equation.  The  last  term  of 

the  parenthesis  becomes  infinite  at  the  stagnation  point  if  is  not 

dx 

here  zero.  In  order  to  avoid  discontinuities  at  the  stagnation  point, 

one  must  therefore  equate  there  =  0.  Then  the  last  term  becomes 

dx 

everywhere  negligible,  and  one  obtains  the  equation 


u 


dc  .  v  dc  ___  a  d2c 
dx  dy  dy2 


which  thus  is  identical  with  the  one  in  the  case  of  two-dimensional  flow. 


1.  Use  of  the  function 


Because  of  symmetry  and  of  the  requirement  ~  =  0  one  uses  for 

,  dx 

p  h 

c  the  expression  c  =  cq  +  c2x  +  c^x  +  .  .  .  with  the  following 
boundary  conditions  for  cy: 


y  =  0;  c0  =  c2  =  c^  =  .  .  .  *  0; 

y  =  °°;  c0  =  c2  =  =  .  .  .  =  0 

For  the  cv  one  obtains  equations  which  contain  rv  and  Uy.  In 

order  to  eliminate  these  constants,  one  may  make  the  following 
substitutions: 


_ _  2U5  /  r3ul  ' 

Ysp*l i  ♦v  -  as  before;  c0  =  F0;  c2  =  —  |G2  +  Bg 


—  Fk  »  —*■  G4  +  \  + 

U1  U1  \  rlu5 


u52  „  ,  r3u3  T  .  r32“l 

-  \  +  -  +  p % 

ulu5  rlu5  rl  u3  / 


Boundary  condition 

rj  =  0;  Fo  =  1;  remaining  functions  =  0; 
tj  *  00;  all  functions  =  0 


IT 


For  the  new  functions  the  following  equations  of  the  second  order 

•  ••  are  obtained: 

t  •  • 

•  •• 

1  ###  —  TT/n  M  —  f*  "P1  1 

»•  ••  :* o  -  ~r2*o 

i  •  •  u 

•  • 

o  g2"  =  -f2G2'  +  f2'G2  “  2g4F0* 

1  H2"  =  -f2H2'  +  f2'H2  -  2h4F0’  +  |  f2F0’ 

i  04"  =  -f204'  +  2f2’G4  -  36gF0' 

1  V  »  +  2f2'H4  -  3h6FQ'  +  }  f/0' 

i  K4"  «  -f2%'  +  2fa'K4  -  5kgF0'  +  |  g4'G2  -  |  g4G2' 

\  V  '  -f2^'  +  2VJ4  *  3W  +  |(VH2  +  VG2)  - 

|(^H2'  +  W)  -  f  (-  f2G2’  +  VG2  -  2g4F0 ') 

|  *  -f2V  *  2f2,94  '  3%Fo'  +  I  h4’H2  -  f  W  - 

|(-f2H2’  +  VH2  -  2h4Fo'  +  I  Vo’ 

2.  Use  of  the  function  >|r 

With  the  same  power  development  for  c  and  the  same  definitions 
of  the  functions  Fq,  G2,  G^,  etc.,  one  obtains 

1  F0”  =  -fjFo’ 

i  g2"  =  -f1G2'  +  fj^’G;,  -  2gjF0' 

5  H2"  .  -flH2'  +  fl'H2  -  2h3F0’  -  i  fjF0- 
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lft."  = 


-flV  +  2f1'G4  -  3g5F0' 


1  TT  " 
“ 

a  ^ 


Ik"  = 


1  J.  " 


i  v- 


-fiV  +  -  5h5Fo'  - 1 W 

-fiV  +  2fiA  -  3W  +  |  e5'Ga  -  f  *>V 

-f^*  +  2f1,J^  -  3J5F0'  +  l^'Hg  +  h^'G^-  I^Gg'  + 
g3Fo’)  "  j(g3V  +  h3G2!) 

-fiV +  a  a  -  Jy0' + 1  yH2  -  f  yy  -  f (qy +  y0')+ 

JfiV 


One  can  show  easily  by  application  of  the  relationships  between  the 
functions  with  even  and  the  functions  with  odd  subscripts  that  the  systems 
of  equations  for  the  cases  1  and  2  are  identical. 

The  first  equations  of  the  two  systems  are  identical  with  the  first 
equation  for  two-dimensional  flow  and  are,  therefore,  also  solved  by 
quadratures . 


C.  FINAL  EXPRESSIONS  FOR  THE  TRANSFER 


The  transfers  are  made  dimensionless  by  the  Nusselt  number 

Nu  =  or,  respectively,  -2-  It  is  easily  shown  that 

Acm  8s8t  Xto  8s8t 

-~=r  =  (-  — |  where  c  and  y  are  "dimensionless  and  without 

v/ite  \  8tj/0  dy 

Reynolds  number."  The  heat  transferred  by  radiation  is,  of  course,  not 
contained  in  this  expression.  For  two-dimensional  symmetrical  bodies  one 
obtains 


=  -  F, 


V/3 


4u,F. 


+  - 


ulu3 
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Corresponding  expressions  are  obtained  in  the  other  cases.  As  one 
can  see  from  the  equations,  one  may  easily  calculate  the  Nusselt  number 
for  arbitrary  pressure  distributions  and  body  shapes  which  agree  with 
the  formulations,  if  one  has  made  a  numerical  calculation  of  the  func¬ 
tions.  Unfortunately,  the  quantity  a  is  left  over  and  one  must  there¬ 
fore  make  different  solutions  for  different  media.  As  is  shown  in  a 
section  below,  however,  one  can  free  the  equations  of  a,  too,  if  a 
is  large. 


D.  NUMERICAL  CALCULATIONS 


For  the  two-dimensional  symmetrical  case  and  for  rotationally 
symmetrical  bodies  the  author  has  numerically  calculated  various  func¬ 
tions,  corresponding  to  the  three  first  terms  of  the  power-series  devel¬ 
opments  in  x.  The  method  of  Runge  and  Kutta  (ref.  14 )  was  used  for 
this  purpose.  This  method  is  rather  time  consuming  but  one  has  good 
possibilities  of  determining  the  errors .  The  first  function  fq  of 
the  two-dimensional  case  has  been  calculated  by  Hiemenz  and  Howarth  with 
an  accuracy  sufficient  for  this  investigation;  Howarth' s  values  are 
directly  used  here.  For  fq  in  the  rotationally  symmetrical  case  there 

exists  a  table  by  Hartree  (ref.  15)  which  was  set  up  by  using  a  mechanical 
differential  analyzer.  The  accuracy  is  here  not  sufficient  and  the  first 
two  derivatives  also  are  required;  for  this  reason  the  function  is  cal¬ 
culated  here  anew.  Since  the  equation  for  fq  is  not  linear  and  one 

therefore  cannot  find  the  solution  by  combining  two  particular  solutions, 
it  was  valuable  to  have  approximate  information  on  f1"  for  t)  =  0. 

The  functions  were  solved  mostly  by  steps  of  q  =  0.2.  Since  the  values 
with  t)- interval  0.1  must  be  known  for  the  successive  calculations,  the 
values  lying  between  were  interpolated  by  means  of  a  Taylor  series.  For 
the  transfer  boundary  layer  the  calculations  for  the  a- value  of  the  air 
(0.7)  were  performed  in  the  two-dimensional  symmetrical  case  because 
experimental  results  for  the  heat  transfer  of  a  circular  cylinder  in  air 
exist  (see,  for  instance,  the  compilation  by  Kroujiline  (ref.  l).  For 
F0")o  one  may  obtain  values  from  a  table  given  by  Goldstein  and  cal¬ 
culated  by  Squire  (ref.  2)  also  in  the  case  of  other  a-values.  Squire 
indicated  an  analogous  expression  for  the  heat  transfer  at  the  stagnation 
point.  For  the  rotationally  symmetrical  case  calculations  have  been 

carried  out  for  a  =  ^  because  the  only  experimental  result  for  the 

transfer  distribution  has  been  found  for  the  evaporation  of  naphthalene 
spheres  (ref.  11),  naphthalene  has  this  value  of  a.  If  one  wants  to 
calculate  the  higher  terms  of  the  power-series  development  in  x  for  a 
special  case,  one  may  combine  the  separate  functions  in  a  single  term  in 
order  to  save  work  expenditure;  but  the  generality  of  the  solution  is  lost 
thereby.  This  has  been  done  for  the  boundary  layer  of  the  body  of  revolution 
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for  a  =  ^  The  parenthesis  in  the  defining  equation  for  G^, 

etc.,  has  teen  combined  into  a  single  function  F^  and  calculated  for 
the  sphere  for  the  pressure  distribution  of  Fage.  (See  below.) 


The  tables  of  the  calculated  functions  axe  printed  at  the  end  of 
the  report  with  the  exception  of  the  higher  ones  for  rotationally  sym¬ 
metric  boundary  layer .  Here  one  has  for  n  -  0 ( and  a  =  — - — ) 

V  0.395/ 

G2'  =  0.3186,  H2'  -  -0.1005,  and  Fl^ '  =  -0.2118.  The  error  of  the 

tabulated  functions  which  will  be  discussed  in  more  detail  later  is  at 
most  a  few  units  in  the  last  digit. 


From  the  tables  one  obtains  for  air,  for  the  pressure  distribution 
measured  by  Hiemens  (ref.  6)  at  Re  ~  19000  for  the  circular  cylinder 


U  -  3- 6314 x  -  2.1709x3  -  1.5l44x5 


-^r  =  0.9449  -  0.5100x2  -  0.5956x\  .  . 

yRi 

The  quantity  x  is  here  dimensionless  (the  length  dimension  x 
divided  by  the  diameter  D) .  Not  only  in  the  range  0°  -  55°  where  the 
series  is  to  apply  exactly  (see  E,  2)  but  up  to  the  separation  point  this 
equation  is  in  good  agreement  with  the  compilation  of  experimental  dis¬ 
tribution  curves  indicated  by  Kroujiline  (ref.  1).  The  derivations  will 
be  discussed  later. 

For  the  sphere  a  qualitative  agreement  with  the  values  obtained  for 
evaporation  of  naphthalene  at  higher  Re  (ref.  11)  is  attained  if  the 
pressure  distribution  according  to  Fage  (ref.  17)  is  used  which  gives 

U  =  3x  -  3.4966x5  +  4.7391*5  -  5 -4181x7  for  Re  =  157200  (ref.  18). 

One  then  obtains 

-7^=  =  1.8615  -  2.l477x2  +  2.4609xi+  .  .  . 

VRe 

The  deviations  depend,  among  other  things,  on  the  fact  that  Fage's 
pressure  distribution  is  possibly  not  fulfilled  for  Reynolds  numbers  as 
small  as  those  used  here.  Later  on  a  more  exact  comparison  will  be  made 
with  more  recent  experimental  values  obtained  by  the  author  at  still 
higher  Reynolds  numbers . 


21 


•  •  • 
•  •  • 
•  •• 

•  • 

•  • 


•  •  •• 

•  •  • 
•  • 


E.  DETERMINATION  OF  THE  ERRORS  OF  THE 
BROKEN-OFF  POWER-SERIES  DEVELOPMENTS 


One  may  use  various  methods:  (l)  The  following  term  is  calculated 
and  for  not-too-large  x  this  indicates  the  error.  (2)  If  the  value 
of  the  errors  is  not  required  with  a  very  high  accuracy,  the  coefficient 
of  the  x-terms  may  be  assumed  to  be  of  the  same  order  of  magnitude.  (3) 
Use  of  a  continuation  method  of  the  profile.  (See  following  section.) 
(4)  In  the  case  of  direct  differentiation,  with  respect  to  y,  of  the 

value  of  —  taken  from  the  transfer  boundary -layer  equation  one  may 
ox 

obtain,  for  transition  to  y  =  0,  the  first  derivative  of  Nu  in  x 
whereby  a  continuation  step  may  be  taken  directly  with  respect  to  Nu. 
Later  on  numerical  calculations  according  to  some  of  these  methods  will 
be  given. 


STEPWISE  DEVELOPMENT  OF  THE  BOUNDARY-LAYER  PROFILE 


A.  TWO-DIMENSIONAL  CASE 
a.  Flow  Boundary  Layer 


Prandtl  (ref.  8)  indicated  for  this  case  a  method  which  is  based 
on  the  fact  that  one  may  obtain  from  the  equations  an  expression  for 


du 


^  containing  only 
ox 


u 


with  derivatives  for  a  prescribed 
with  dimensionless  quantities  without  Reynolds  numbers 


~  becomes 
ox 


du 

dx 


+ 


When  two  adjacent  profiles  (at  x  -  Ax  and  x)  were  known,  for 
instance,  by  application  of  the  method  of  Blasius  and  Hiemenz,  it  was 
possible  to  calculate  a  third  for  x  *•  Ax.  With  the  u-values  at  x  -  Ax 

the  2Ax  r—  values  for  x  were  used.  One  could  then  continue  in  the 
dx 

same  manner  with  the  profiles  at  x  and  x  +  Ax.  In  order  to  guarantee 
the  convergence  of  the  expression,  one  was  not  to  use  the  original  numer¬ 
ical  profile  at  x  but  had  to  replace  it  by  another  which  satisfied 
certain  requirements.  In  order  to  calculate  those,  u  was  developed 
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\ —  ^vyv 

into  a  power  series  with  respect  to  y:  u  =  — r~*  For  the  ®v  one 

v=l  v* 

obtained  certain  conditions  by  substitution  into  the  flow  equation  whereby 
only  some  of  them  could  be  chosen  arbitrarily.  These  latter  were  deter¬ 
mined  by  comparison  with  the  given  profile.  Gortler  (ref.  9)  perfected 
the  method  practically  and  used  it  in  Hiemenz'  pressure  distribution  over 
the  circular  cylinder.  In  the  present  report  corresponding  ideas  are 
used  for  bodies  of  revolution  and  for  the  transfer  boundary  layer,  and 
the  necessary  expressions  are  added  and  discussed. 


b.  Transfer  Boundary  Layer 
From  the  basic  equation  one  obtains  directly 

1/1  d2c 


dc 

bx 


—  v 

u\o  by2  by) 


This  equation  may  therefore  be  used  directly  for  step-by-step 
continuation  of  the  vapor  and  temperature  boundary  layer.  Conditions 
become  here  simpler  insofar  as  no  integration  is  necessary.  However, 
here  also  the  danger  exists  that  the  expression  becomes  uncertain  at  the 

be 


wall  (because  of  u  occurring  in  the  denominator).  Moreover, 


dx 


must 


become  identically  zero  at  the  wall.  In  order  to  satisfy  the  require¬ 
ments,  one  resolves  here  also  the  quantity  c  into  a  power-series  devel¬ 
opment  with  respect  to  y(bv  function  of  x  only): 


c  =  1 


-Z 


fcvyv 


v=l  v* 


By  substitution  one  obtains 


ay,yV'Hl  y~WyV~1 
(v  +  (v  -  i): 


ly — bvyv~2 

°^_(v  -  2)! 


By  comparison  of  terms  of  the  same  degree,  one  arrives  at  the 
relation  between  the  bv.  For  the  first  nine  bv  there  applies  (with 

f  =  -UU') 
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bp  free;  b2  =  =  0 

=  2a1b1'  -  a1'b1  free;  =  3^'  -  f  'b-^  b^  =  0 
^  =  10a12ab1"  +  5a]a1,(l  -  3a)b1'  +  ja^^lOa  -  l)  -  W'^a  +  1))^  free 

^  =  48a1ofb1"  +  2 J2a_^f '  ( 3  -  7a)  -  l^ap'fajb-L'  +  [-2aif"(3a  +  1)-  15fai"a  + 

f'a1l(35a  -  2)] bl5 

^  =  63f2ab1"  +  7ff'(2  -  9a)b1'  +  [f'2(35a  -  2)  -  ff"(2lcr  +  2)]  bp 


The  free  coefficients  are  calculated  as  before  by  comparison  with 
the  given  profile,  and  the  c -values  developed  in  power  series  are  sub- 

stituted  into  the  above  equation  for  — . 

ox 


B.  ROTATIONALLY  SYMMETRICAL  CASE 


a.  Flow  Boundary  Layer 


Hie  equations  read 

u^u+v^u=UU. 

ox  Oy  byd 

d(ur)  +  d(vr)  =  Q 
bx  by 

Here  r  signifies,  as  before,  for  a  blunt  body  of  revolution  the 
distance  between  the  axis  of  rotation  and  the  base  point  of  the  normal 
to  the  surface. 

By  eliminating  ~  one  obtains  a  linear  equation  of  the  first  order 
ox 

in  v  with  solution 


r  UU’  1  bTv\  ,  1  dr 

v  =  -u  /  -  + - dy - - —  uy 

Jo  L2  u2  3y2/  r 
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By  forming  the  derivative  from  v  with  respect  to  y  and  using 

,•/*,  the  continuity  equation  one  obtains 

•  •• 

•  •  _ 

+  1  d£.  y  bu 
r  dx  by 


For  given  u-profile,  one  may  therefore  use  an  equation  for  the 
continuation  of  the  boundary -layer  profile  which  differs  from  the  equa¬ 
tion  for  the  two-dimensional  case  only  with  respect  to  the  last  term. 

In  order  to  establish  the  convergence  at  the  wall,  here  also  a  power- 
series  development  in  y  becomes  necessary 


•  •• 
•  •  •  • 


du  _  d 
dx  by 


b£ u 
'by2j 


|hyj 


since 


u 


£H,yv 

~VT 


>\  „v+l 


g 


V  .  -  £L-  ♦  a„  - 

•  Mv  r  J(v  +  i): 


By  substitution  into  the  basic  equations  one  obtains  ^with  f  = 


-UU'; 


free;  a2  =  f;  a^  =  0; 

a^  =  aipi'  -  a]_2g  free;  as  =  2a1f'  -  4a-|_fg;  a^  =  2ff'  -  4f2g 


As  before,  one  determines  the  free  coefficients. 

b.  Transfer  Boundary  Layer 

With  the  same  expression  as  above  for  the  two-dimensional  case  one 

obtains  for  the  bv 


oc. 


•  •  • 

•  •• 

•  • 

•  • 

•  ••• 

•  •  •• 

•  •  • 

•  • 


b-^  free;  =  b^  -  0; 

b)  b^ 

=  Za^'  -  a1'b1  -  a1b1g  free;  =  3?\'  -  f'\  -  ?\&i  \  '  °? 


The  practical  execution  in  the  last  three  cases  will  be  discussed 
in  a  later  report  together  with  numerical  calculations.  The  methods  of 
continuation  discussed  yield  results  the  accuracy  of  which  depends 
exclusively  on  the  work  expenditure  and  is  therefore  not  limited  by 
postulating  approximation  functions.  The  methods  may  also  be  used  for 
determining  the  accuracy  of  the  aforementioned  power-series  developments 
in  x  in  the  case  of  breaking-off  after  a  certain  number  of  terms  at  a 
certain  point.  One  then  starts  the  continuation  method  at  an  x  so 
small  that  the  error  is  certainly  small,  and  compares  the  result  then 
obtained  at  a  larger  x  with  the  one  directly  calculated  from  the  power- 
series  development  in  x. 


DEPENDENCE  OF  THE  EVAPORATION  AND  THE  HEAT  TRANSFER  ON  a 

A.  GENERALITIES 


Pohlhausen  (ref.  10)  has  shown  for  the  plane  that  Nu  is  approxi¬ 
mately  proportional  to  the  quantity  ^/~a.  In  the  approximate  calcula¬ 
tions  of  Kroujiline  (ref.  1)  the  same  was  shown  for  the  circular  cylinder. 
Ulsamer  (ref.  10)  demonstrated  that  the  law  may  be  approximately  selected 
from  various  experimental  investigations  on  the  heat  transfer  of  a 
circular  cylinder.  The  author  of  this  report  has  confirmed  the  law  at 
least  approximately  in  the  case  of  evaporation  of  drops  (ref.  11). 

From  the  equations  of  the  section  on  power-series  developments  in 
x  one  sees  that  a  can  probably  not  be  eliminated  from  them  by  simple 
transformations.  Thus  one  cannot  expect  a  relation  as  simple  as  the 
aforementioned  to  apply  exactly.  For  the  case  where  a  is  very  large 
and  the  transfer  boundary  layer  therefore  thin  compared  to  the  flow 

boundary  layer,  the  author  of  this  report  found  the  \J  o-law  to  be  exact. 
In  this  case  the  curvature  of  the  velocity  profile  may  be  neglected  in 
the  entire  transfer  boundary  layer,  and  one  may  replace  u  by  (u1)^ 

and  v  by  ^(v")Qy^  in  the  general  boundary -layer  equation,  with  the 
apostrophes  indicating  derivatives  with  respect  to  y. 
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(u.)c,|t<z2oy2|  =  ig 

The  variable  £  =  y  ^ fa  is  introduced  for  c  (not  for  u  and  v' 


/  ,  %  6c  ^0  »2  6c  62c 

(u)rf£  +  — 5  ac’p 


Boundary  conditions:  £  =  0;  c  =  1.  £  =  c  =  0. 

Thus  one  has  obtained  an  equation  free  of  a.  For  this  reason,  c 


becomes  c  =  f[x,y^J~oj;  hence  follows  that  for  large  a  the  quantity 

Nu  is  proportional  to  the  quantity  \/o  on  the  entire  surface  in  the 
boundary  layer.  That  the  same  law  has  been  found  experimentally  also 
for  a  a  that  is  not  large,  is  based  on  the  fact  that  the  quantity 

Nu  If  does  not  vary  greatly  with  a  and  may  therefore  be  found  to  be 
approximately  constant  in  a  small  region. 


B.  TWO-DIMENSIONAL  SYMMETRICAL  CASE 


Into  the  equations  for  Fq,  Fg,  G4,  %  .  .  .  the  following 

functions  and  variables  are  introduced: 

5  =  1  \/f(fi")o*  Vn)  -  r2<i)  =  MU 

Vb  (fl  )o 


(g5")o  (V)n 

G4  =  irr  r4(0;  H4  =  j~-nv- -  04(1);  etc 

( 1 )°  (L  )o 


Boundary  conditions: 


i  =  0;  <I>q  =  lj  remaining  functions  =  0; 
|=°°;  all  functions  --  0. 
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Taking  into  account  that  for  large 
is  valid,  one  obtains 


a  the  equation  \|r 


!<*")oy2 


V 


ft 


-3t2V 

-34  2  $2 ' 

+  12 ^ $2 

-342r4' 

+  24^r4 

-342ei+' 

+  24404 

94 2 V 
15l^0' 

154 2$0'  + 


For  the  pressure  distribution  according  to  Hiemenz  (see  p.  20) 

Nu 

one  obtains  in  the  case  of  a  circular  cylinder  -—■r^.-==.  =  1.2592  - 

VO\pe 

0.7583x2.  .  .  ;  for  the  case  calculated  above  c  =  0.7  one  obtains  for 

the  corresponding  quantity  1.0642  -  0.5744x2  -  0 . 6708x1+  ....  Here 
x  signifies  the  dimensionless  abscissa  which  is  obtained  from  the  length 
dimension  through  division  by  the  diameter  D.  The  functions  Oq  and 
$2  given  numerically  in  table  6. 


C.  ROTATIONALLY  SYMMETRICAL  CASE 


The  following  functions  are  now  introduced: 


^  ^  \[^  1  )o5  F0^  “  G2^) 


e2^) 


The  equations  become 


V  =  -36^o' 

r2"  =  -342r2'  +  6\ r2  -  6|^0* 

e2"  .  -}(se2'  +  6ie2  -  6^0'  -  §  52*Q 

(h3  )o 
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As  in  the  previous  case,  the  solution  of  the  first  equation  is 


^  =  1 


JO 


For  the  pressure  distribution  according  to  Fage  (ref.  17)  (see  p.  20) 


one  obtains  for  the  sphere 
1 


Nu 


\f4 


=  =  1.4723  -  ( 


•)x‘ 


for 


Re 


a  = 


0.395 


one  obtains  1.3658 


DISCUSSION  OF  APPROXIMATE  METHODS 


As  has  been  mentioned  above,  Pohlhausen  (ref.  3)  gave  an  approxi¬ 
mate  method  for  the  solution  of  the  boundary -layer  equation  for  the  flow 
about  a  circular  cylinder.  Tomotika  (ref.  18)  applied  this  method  to 
the  sphere.  Kroujouline  used  a  similar  method  for  the  transfer  for  the 
cylinder,  applying  a  broken-off  series  development  in  y  which  was 
determined  with  utilization  of  the  integral  condition  of  the  transfer 
boundary  layer.  For  the  flow  boundary  layer  he  used  a  parabolic  profile 
whereby  the  agreement  may  be  assumed  to  be  bad  particularly  in  the  case 
of  pressure  increase.  Probably  better  approximations  could  have  been 
obtained  with  the  use  of  polynominals  of  the  fourth  degree.  These 
statements  are  valid  only  when  the  transfer  boundary  layer  is  thinner 
than  the  flow  boundary  layer.  Here  a  brief  description  is  given 
concerning  some  considerations  of  the  author  of  this  report  concerning 
a  body  of  revolution,  for  various  relative  magnitudes  of  the  two  layers. 

The  integral  condition  formerly  not  set  up  for  bodies  of  revolution 
becomes  (see  p.  4) 


1  d_ 
r  dx 


which  may  be  derived,  for  instance,  by  integration  of  the  original 
equation.  Here  5  is  the  thickness  of  the  transition  boundary  layer. 
Using  dimensionless  quantities  without  Reynolds  numbers  only,  A  is 

replaced  by  This  is  assumed  below. 


If  for  the  two  profiles  the  definitions 


and 


c  ,  !  -  2  *+  2!*V  -  (t 

6  \b)  U 


are  used,  one  obtains  for  5  <  6]_  the  equation 


2r  _  d_ 
5nZ0  dx 


/a-,z2  a?z^  3a,  z^  ai,z^i 
rUbn  \-±—  +  -f —  +  —3 —  +  -=t_— J 
1 V 15  42  280  180/ 


Here  6^  is  the  thickness  of  the  flow  boundary  layer  and  z  =  . 


For  6  >  6p,  the  integration  is  performed,  with  use  of  the  integral 
condition,  first  from  0  to  6^  and  then  from  6-^  to  6.  Result 


2r 

5-^zo 


d_ 

dx 


rUh 


3£  +  (!l  +  !£  ♦  !i  +  .  j]  .  2;^  +  f2  +  !i  +  \  . 


10  '  2  ’  3  ’  4  ’  5  7  z '  3  ‘  4  '  5  ’  6 

2  i  +  +  h  +  .  i'l  .  Lb.  +  tz  +  S  +  !t .  i] 


8 


r4  \6 


8 


) 


The  two  equations  have  the  same  form  when  the  parentheses  after 
6l  are  denoted,  for  instance,  by  the  letter  P.  From  the  first  of  the 

two  equations  one  sees  that,  for  a  a  so  large  and  a  z  therefore  so 
small  that  only  the  first  term  of  the  parenthesis  must  be  considered, 
this  z  is,  for  a  given  x, 


k,  inversely  proportional  to  the  quantity 


Since  the  Nusselt  number  Nu  equals 


2 —  \/Re, 
)]_Z  V  ’ 


Nu  is,  for  a 


large  a,  proportional  to  <?fa  also  according  to  this  approximate 
theory .  V 


Since  r,  U,  a v,  and  6^  are  known  functions  of  x,  we  have  in 
any  case  an  equation  of  the  first  order  with  z  and  x  which  can  be 
solved  with  customary  methods  (for  instance,  with  the  isocline  method 
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or  according  to  Runge  and  Kutta).  The  only  boundary  condition  required 
for  this  is  the  z-value  at  x  =  0.  This  value  is  calculated  from  the 
equation  AzcP  =  1  where  A  =  U'5^2  is  identical  with  the  parameter  A 
used  by  Pohlhausen  and  Tomotika.  For  the  sphere  where  U’  -  3  and 
A  -  4.716,  one  obtains  in  the  proximity  of  the  stagnation  point 


-  -  O.8759Z5  -  0.2648z^  +  0.0l809z5  +  0.0056lz6 


or,  respectively, 


1  =  1.4l48z2 
a 


1.2295z  +  0.5052 


0.0746  O.OI89 

z  ^ 


For  a  given  z  and  therefore  also  given  Nu 
calculate  the  corresponding  a.  For 


one  may  easily 


z  =  0.0 

4  $  - 0 


0.1  0.4  0.7  1.1  1.6  2.0  3-0  4.0,  one  obtains 

1.526  1.511  1.464  1.4l8  1.356  1.284  1.232  1.128 

1.049  and 

0.095  0.367  0.622  0.935  1.288  1.545  2.121  2.631 


The  quantity  Nu j  \^Re  \fa  is  therefore,  for  a  large  a,  almost 
constant  and  varies  in  the  proximity  of  the  stagnation  point  about 


linearly  with 


The  reason  for  choosing,  above,  z  instead  of  .6  as  the  dependent 
variable  was  that  z  probably  varies  little  with  x  (compare  Kroujouline 
(ref.  l))  and  can  therefore  be  calculated  exactly  more  easily. 


In  the  later  more  detailed  report  on  the  investigations,  the  numer¬ 
ical  results  of  this  formulation  as  well  as  of  others  will  be  discussed. 
It  was  shown  that  the  choice  of  the  profile  form  had  a  great  effect  on 
the  result. 


SUMMARY 


A  preliminary  report  is  given  of  a  theoretical  investigation  of  the 
boundary -layer  flow  for  two-dimensional  and  rotationally  symmetrical 
bodies.  The  evaporation,  the  heat  transfer,  and  the  velocity  are  cal¬ 
culated  by  power-series  developments  with  respect  to  the  meridian  length. 


The  coefficient  functions  which  were  calculated  numerically  in  some  cases 
have  been  chosen  so  that  the  calculation  is  valid  for  all  pressure  distri¬ 
butions  and  body  shapes.  The  methods  for  determination  of  the  errors  in 
breaking  off  the  series  are  briefly  treated.  Methods  of  continuation  are 
discussed.  It  is  shown,  for  large  Prandtl  numbers,  that  the  Nusselt 
number  is  exactly  proportional  to  the  cube  root  of  the  Prandtl  number. 
Finally,  approximate  methods  of  calculation  are  discussed. 
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TABLE  2.  TWO-DIMENSIONAL  SYMMETRICAL 
TRANSITION  BOUNDARY  LAYER 


n 

1  -  F0 

-Ftf' 

F2 

f2’ 

G4 

G4  ' 

H4 

V 

0.0 

0 

0.4959 

0 

-0.1119 

0 

-0.0977 

0 

0.0318 

.2 

.0991 

.4953 

-.0224 

-.1113 

-.0195 

-.0970 

.0064 

.0320 

.4 

.1979 

.4917 

-.0443 

-.1077 

-.0386 

-.0932 

.0129 

.0338 

.6 

.2954 

.4825 

-.0651 

-.0988 

-.0565 

-.0846 

.0200 

.0378 

.8 

.3904 

.4660 

-.0835 

-.0841 

-.0721 

-.0704 

.0281 

.0433 

1.0 

.4813 

.4416 

-.0983 

-.0638 

-.0843 

-.0517 

.0374 

.0490 

1.2 

.5 666 

.4095 

-.1087 

-.0397 

-.0926 

-.0301 

.0476 

.0529 

1.4 

.6447 

.3708 

-.ll4l 

-.0138 

-.0964 

-.0078 

.058’ 

.0530 

1.6 

.7146 

•  3275 

-.1143 

.0111 

-.0958 

.0131 

.0685 

.0480 

1.8 

.7755 

.2818 

-.1099 

.0328 

-.0913 

.0307 

.0771 

.0376 

2.0 

.8273 

.2360 

-.1015 

.0496 

-.0838 

.0440 

.0832 

.0228 

2.2 

.8701 

.1924 

-.0904 

.0605 

-.0741 

.0521 

.0861 

.0054 

2.4 

.9045 

.1526 

-.0777 

.0656 

-.0633 

.0554 

.0854 

-.0123 

2.6 

.9315 

-1177 

-.0645 

.0654 

-.0522 

.0544 

.0813 

-.0281 

2.8 

.9520 

.0883 

-.0518 

.0610 

-.0417 

.0502 

.0744 

-.0404 

3-0 

.9672 

.0644 

-.0403 

.0540 

-.0323 

.0440 

.0654 

-.0481 

3.2 

.9781 

.0457 

-.0303 

.0454 

-.0242 

.0368 

.0555 

-.0511 

3.4 

.9858 

.0315 

-.0221 

.0366 

-.0176 

.0295 

.0453 

-.0500 

3.6 

.9910 

.0211 

-.0156 

.0283 

-.0124 

.0227 

.0356 

-.0458 

3.8 

.9944 

.0138 

-.0107 

.0211 

-.0084 

.0168 

.0271 

-.0396 

4.0 

•  9966 

.0088 

-.0071 

.0151 

-.0056 

.0119 

.0199 

-.0325 

4.2 

.9980 

.0054 

-.0046 

.0104 

-.0036 

.0082 

.0141 

-.0254 

4.4 

.9989 

.0032 

-.0029 

.0070 

-.0022 

.0055 

.0097 

-.0191 

4.6 

.9994 

.0019 

-.0018 

.0045 

-.0014 

.0035 

.0064 

-.0137 

4.8 

.9996 

.0011 

-.0010 

.0028 

-.0008 

.0021 

.0041 

-.0095 

5.0 

•  9998 

.0006 

-.0006 

.0017 

-.0005 

.0013 

.0026 

-.0063 

5.2 

.9999 

.0003 

-.0003 

.0010 

-.0003 

.0008 

.0015 

-.0040 

5.4 

1.0000 

.0002 

-.0002 

.0006 

-.0002 

.0005 

.0009 

-.0025 

5-6 

.0001 

-.0001 

.0003 

-.0001 

.0003 

.0005 

-.0015 

5.8 

.0000 

-.0000 

.0002 

-.0000 

.0002 

.0003 

-.0009 

6.0 

.0001 

.0000 

.0001 

-.0005 

6.2 

.0000 

.0001 

-.0002 

6.4 

.0000 

-.0001 

6.6 

-.0001 

6.8 

-.0000 

TABLE  3.  ROTATIONALLY  SYMMETRICAL  FLOW 
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BOUNDARY  LAYER 


B 

fl 

*l’ 

V 

— 

«3 

83' 

83" 

h3 

h3' 

h5" 

0.0 

0 

0 

0.9277 

0 

0 

1.0475 

0 

0 

0.0448 

.1 

.0046 

.0903 

.8777 

.0051 

.0998 

.9477 

.0002 

.0044 

.0448 

.2 

.0179 

.1755 

.8277 

.0196 

.1896 

.8488 

.0009 

.0090 

.0444 

.3 

.0395 

.2558 

.7778 

.0427 

.2696 

•  7517 

.0020 

.0133 

.0434 

.4 

.0689 

.3311 

.7282 

.0732 

.3400 

.6574 

.0036 

.0176 

.04 16 

.5 

.1056 

.4014 

.6788 

.1104 

.4012 

.5666 

.0055 

.0217 

.0391 

.6 

.1490 

.4669 

.6300 

.1532 

.4535 

.4802 

.0079 

.0254 

.0356 

•  7 

.1988 

.5275 

.5319 

.2008 

.4974 

.3986 

.0106 

.0288 

.0314 

.8 

.2544 

.5833 

.5348 

.2524 

.5334 

.3227 

.0136 

.0316 

.0265 

.9 

•  3153 

.6345 

.4888 

.3072 

.5621 

.2528 

.0169 

.0340 

.0210 

1.0 

.3811 

.6811 

.4443 

.3646 

.5842 

.1895 

.0204 

.0358 

.0152 

1.1 

.4514 

.7234 

.4014 

.4239 

.6002 

.1328 

.0241 

.0370 

.0091 

1.2 

.5256 

.7614 

.3604 

.4845 

.6110 

.0832 

.0278 

.0377 

.0032 

1.3 

.6035 

•  7954 

.3215 

.54  59 

.6171 

.0403 

.0316 

.0377 

-.0026 

1.4 

.6846 

.8258 

.2850 

.6078 

.6193 

.0044 

.0353 

.0372 

-.0080 

1.5 

.7686 

.8526 

.2508 

.6696 

.6182 

-.0251 

.0390 

.0361 

-.0127 

1.6 

-8550 

.8761 

.2192 

.7313 

.6144 

-.0483 

.0425 

.0346 

-.0168 

1.7 

•  9437 

.8966 

.1901 

.7925 

.6087 

-.0657 

.0459 

.0327 

-.0202 

1.8 

1.0342 

.9142 

.1637 

.8530 

.6015 

-.0780 

.0491 

.0306 

-.0228 

1.9 

1.1264 

.9294 

.1398 

.9127 

.5932 

-.0857 

.0520 

.0282 

-.0244 

2.0 

1.2200 

.9422 

.1185 

•  9716 

.5845 

-.0894 

.0547 

.0258 

-.0254 

2.1 

1.3148 

•  9530 

.0996 

1.0296 

.5755 

-.0898 

.0572 

.0233 

-.0256 

2.2 

1.4106 

.9622 

.0831 

1.0867 

.5 666 

-.0876 

.0594 

.0207 

-.0252 

2.3 

1.5072 

.9698 

.0688 

1.1430 

.5580 

-.0834 

.0613 

.0182 

-.0243 

2.4 

1.6045 

.9760 

.0564 

1.1984 

.5500 

-.0776 

.0630 

.0158 

-.0229 

2.3 

1.7024 

.9811 

.0458 

1.2530 

.5425 

-.0709 

.0645 

.0136 

-.0212 

2.6 

1.8007 

.9853 

.0370 

1.3069 

.5358 

-.0637 

.0657 

.0116 

-.0193 

2.7 

1.8994 

.9886 

.0296 

1.3602 

.5298 

-.0563 

.0668 

.0097 

-.0174 

2.8 

1.9984 

.9912 

.0234 

1.4129 

.5245 

-.0490 

.0677 

.0081 

-.0153 

2.9 

2.0977 

.9932 

.0184 

1.4651 

.5200 

-.0420 

.0684 

.0067 

-.0133 

3.0 

2.1971 

.9949 

.0143 

1.5169 

.5161 

-.0356 

.0690 

.0054 

-.0114 

3.1 

2.2966 

.9962 

.0110 

1.5683 

.5128 

-.0297 

.0695 

.0044 

-.0097 

3.2 

2.3963 

.9972 

.0085 

1.6195 

.5102 

-.0245 

.0699 

.0035 

-.0082 

3.3 

2.4961 

•  9979 

.0064 

1.6704 

.5079 

-.0200 

.0702 

.0028 

-.OO67 

3.4 

2.5959 

.9985 

.0048 

1.7211 

.5061 

-.0161 

.0705 

.0022 

-.0054 

3.5 

2.6958 

.9989 

.0036 

1.7716 

.5047 

-.0128 

.0706 

.0016 

-.0044 

3.6 

2.7957 

.9992 

.0026 

1.8220 

.5036 

-.0101 

.0708 

.0013 

-.0035 

3.7 

2.8956 

•  9995 

.0020 

1.8723 

.5027 

-.0078 

.0709 

.0010 

-  .0028 

3.8 

2.9956 

.9996 

.0014 

1.9226 

.5020 

-.0060 

.0710 

.0007 

-.0021 

3.9 

3.0955 

•  9997 

.0010 

1.9727 

.5015 

-.0046 

.0711 

.0006 

-.0016 

4.0 

3.1955 

.9998 

.0007 

2.0229 

.5011 

-.0034 

.0711 

.0004 

-.0013 

4.1 

3.2955 

•  9999 

.0005 

2.0730 

.5008 

-.0026 

.0711 

.0002 

-.0009 

4.2 

3.3955 

•  9999 

.0004 

2.1230 

.5006 

-.0019 

.0712 

.0002 

-.0007 

4.3 

3.4955 

.9999 

.0003 

2.1731 

.5004 

-.0014 

.0712 

.0002 

-.0006 

4.4 

3.5954 

.9999 

.0002 

2.2231 

.5003 

-.0010 

.0712 

.0002 

-.0004 

4.5 

3.6954 

1.0000 

.0001 

2.2731 

.5002 

-.0007 

.0000 

-.0002 

4.6 

.0001 

2.3231 

.5001 

-.0005 

-.0002 

4.7 

.0000 

2.5732 

.5001 

-.0004 

-.0001 

4.8 

.5000 

-.0002 

-.0001 

4.9 

-.0001 

-.0000 

5.0 

-.0001 

5.1 

-.0001 

5.2 

, 

_ 

-.0000 

FLOW  BOUNDARY 
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TABLE  5.  ROTATIONALLY  SYMMETRICAL 
TRANSFER  BOUNDARY  LAYER 


a 

o 

t 

0.5 

0.4129 

•  7 

.4705 

1 

•  5390 

1/0.395 

•  7599 

10 

1.2389 

100 

2.7365 

TABLE  6.  TWO-DIMENSIONAL  SYMMETRICAL  TRANSFER 
BOUNDARY  LAYER,  a  LARGE 


1 

1  -  <&o 

-ft)' 

* 2 

42' 

0.0 

0 

1.1198 

0 

-0.4799 

.1 

.1120 

1.1187 

-.0479 

-.4780 

.2 

.2235 

1.1109 

-.0952 

-.4647 

.3 

-3337 

1.0900 

-.1401 

-.4293 

.4 

.4409 

1.0504 

-.1801 

-.3637 

•  5 

.5430 

.9883 

-.2118 

-.2647 

.6 

.6378 

.9023 

-.2320 

-.1361 

•  7 

.7228 

.7947 

-.2384 

.0099 

.8 

.7962 

.6711 

-.2301 

.1541 

•  9 

.8567 

.5402 

-.2083 

.2747 

1.0 

.9043 

.4120 

-.1765 

.3530 

1.1 

•  9396 

.2959 

-.1395 

.3794 

1.2 

.9641 

.1989 

-.1023 

.3565 

1.3 

.9801 

.1245 

-.0693 

.2981 

1.4 

.9897 

.0720 

-.0432 

.2232 

1.5 

.9951 

.0383 

-.0246 

.1500 

1.6 

.9979 

.0186 

-.0128 

.0904 

1-7 

.9992 

.0082 

-.0060 

.0488 

1.8 

•  9997 

.0033 

-.0026 

.0236 

1-9 

•  9999 

.0012 

-.0010 

.0102 

2.0 

2.1 

2.2 

2.3 

2.4 

1.0000 

.0004 

.0001 

.0000 

-.0003 

-.0001 

-.0000 

.0040 

.0014 

.0004 

.0001 

.0000 

